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Note on the Central Differential Equation in the Relativity 

Theory of Gravitation, 

By Prof. A. E. Forsyth, F.E.S., Imperial College, South Kensington. 

(Eeceived March 4, 1920.) 
The critical equation in Prof. Einstein's relativity theory of gravitation is 



given in the form*' 






where m is the mass of the sun, % is the reciprocal of the distance of the 
moving body from the sun, 

a is the " major semi-axis of the orbit," the units being astronomical, and e is 
the " excentricity " of the orbit. 

The foregoing equation can be integrated exactly, the functions involved 
being elliptic and not circular ; and the approximations for the solar system 
are easily obtained, because the modulus of the elliptic functions is small. 
The necessary analysis is as follows. 

1. The equation is 

\d^) - ah?^h^'' u~\-lmu 

= A^.H3B^62_^3C^6 + I) = U, 
where A = 2m, B=-i, C = ?~ D= ^ 



3 h^' ah?' 

The equation U = has three roots ; and their character, as to reality, 
depends upon the sign of G^+4H^, where (with the usual notation) 

H = AC— B2= -i-f^ — 

9^3 A2' 

G- = A2D-3ABC + 2Bs= -_A-f ^ !!!l-4 ''^ 



27 3 Ih' aK^' 
On substitution and reduction, we have 

C^-l-4H^ = — ^— — m^ 32 m^ I6m^ 256 m^ 

" li^\ 27'F"^27'^ y^I^ '^~h?'^'W^ 

■^ Prof. Eddington's <Eeport on the Eelativity Theory of Gravitation,' 2nd edition, 
p. 50 (1920). 
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By using the relation 

the last expression leads to the equation 

In the solar system, as regards even the nearest known planets, mPjh^ is of 
the order 10"^^, while e^ is never of so small an order ; consequently, the first 
term within the bracket governs the sign of G^ + 4:ffi Thus G^ + 4H^ is 
negative, but not attaining a limit zero within the system ; therefore, by the 
usual rule, the three roots of U = are real, and no two of them are equal. 
Let them be denoted by a, ^, 7, arranged in descending order of magnitude ; 
their product is l/2ah^, a positive quantity ; hence either all three real roots 
must be positive, or one of them must be positive and the other two must be 
negative. 

2. As a matter of fact, the three roots are positive for the solar system. 
Let Aw + B = t ; 

then t^ -\' SKt -\- G = 0. 

Comparing this equation, as usual when the three roots are real, with the 

equation 

cos^^— f cos^— |cos30 = 0, 

we find t =: 2(-— B[)*cos^, 

where cos 30 = ~«.i ; — tt^tt, 

giving a positive value for cos 36, We shall assume a positive value for 30, 
between and |-7r; there is no loss of generality in the assumption. The 
three values of t correspond to the values 0, 0-i-|7r, B-^^tt; hence, as 

cos (04- f tt) = — I cos 0—^^ sin 6, 

/3 
cos (O + ^w) ^ — I cos 6 + ^^r-~ sin 0, 

2 



we have 



Again, 



2mci—^ = 2( — H)^cos0, 
2m/S--^ = 2(-H)icos(0 + f tt), 
2m7-:-i == 2(--H)»eos(0+|7r). 



sin^ 30= 1+ --rppK = 108 — — " 

4W h^ 
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approximately, for the solar system ; so that 

approximately. Also, approximately, 

and therefore 

2my3-i=: 2(^-25') (-K^^) 

1 _1_ 9 _L, — -— 

approximately. Thus ^ is positive. Consequently 7 is positive, and its 
value is given by 

2 W? fV — • -i ^^Z M l i n- JL -t- 2 • =i 



».l 



. Eeturning now to the original differential equation, we have 

Now du/dif} is real For the solar system, 2ma is approximately equal to 1, 
while m = 1-47, and u is of order not higher than 10"^. Hence 

^ being the value of % at perihelion. We shall assume that ^ is 
measured positively throughout the orbit, so that f?^ is positive. From 
perihelion to aphelion %v decreases, so that d% is negative in this range. 
Thus du/d(f> is negative, and therefore . 

Substitute 

u = 7 -f O— 7) sin^ ifr ; 

after a slight reduction, we find 



d^ rsz — 2p 



d'^ 



(1 — ^sin^*^)* 



where p = {2m(a— 7)}"^, 7^2 _.P — 7 

« — 7 

The value of '?// at perihelion is /8, and it subsequently decreases • hence the 
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value of i|r at perihelion is |-7r, and it subsequently decreases. Consequently, 
at any place before aphelion in the orbit, we have 

^P ~~^ J vV (1 —^^^ sin^ l/r)i 



\ 



^ (1 — Psin^i/r)2 



JTT 



^'^ f^ <^^ 



jo(l'-/u^si] 



( 1 — F sin^ '\{r)^ J ( 1 -— /^^ sin^ '^)t 

where <^ = w at perihelion, and, in the usual notation of elliptic functions, 
K is the complete elliptic integral of the first kind, and 

that is, y^r = am /K — ^--^j. 

Now u = 7-f (/3-~ 7).siix^i|r 

. fQ s. l — cn 2ll 
= 7 + (/3— 7) 



^ — w 



= 7 + (/S"~7) 



\ p 



1 4- ^?'??' 



Iz:?- 



which is the accurate equation of the orbit represented by the original 
differential equation. It involves elliptic functions, instead of the circulai: 
functions that occur in the equation of the ISTewtonian orbit. 

4. Further, the value of i|r (decreasing from l^ '^ ^^ perihelion) is zero at 
aphelion ; so, if ^^ denote the value of ^ at aphelion, we have 

From aphelion to perihelion we have similar relations ; so that, if at the 
next perihelion, <^2 be the value of <j[>, which is to begin with <^i at aphelion, 
we have 

^2 — ^1 = 2/?K. 

Thus ^2— tar = 4/3K:. 
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Let the advance of the perihelion in one revolution be denoted by A ; then 

A = {(j}2-—m) — 27r = 4pK— • 27r. 

The expression is accurate, and not merely approximate, in connection with 
the original differential equation so long as the magnitudes in that equation 
allow three real roots for U = 0, whatever be the gravitational system. 

The preceding results are accurate, and not merely approximate, subject to 
the condition just stated. In the exposition of the application to the solar 
system, and specially with reference to the advance of the perihelion of 
Mercury in one complete revolution, approximation is made with regard to 
the magnitudes of the solar system that are involved. To bring the 
approximation into relation with the foregoing accurate results, we must deal 
with the two relations 

^ = ry-f (/S — ry)-- ,^^, A = 4|r)K ~ 27r. 

p 

5. As regards the approximation to A, we require the most important 
terms in p and in K, 
The value of p is 

From the foregoing values of a, yS, 7, we have 

2m(a'-y) = 2(-H)i{cosi9— co8(^ + |7r)} 

= 2(--H)i(|cos^ + -^^sin^' 



3 A h?l2\ ' v^3. 

~-^^^3 h'' 

to the order of approximation that includes m^jh^ but no higher power of 
mill. Hence, to this order, 

P 



The value of F is 



— 


l.sg- 


1 




/8-7 





a — « 



Again from the foregoing values of /3 and 7, we have 
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to the order retained ; and therefore, to this order, 

"ill 



But K = 






*hn 



(1— Psin^'\/r)5 



Jo 

to the order of approximation : that is, to the order, 

K = |-7r(l-f|-F) 

= 27r(l + ~^). 
Thus A = 4/3K-27r 

and therefore there follows Prof. Einstein's inference from his theory that, to 
the order of approximation adopted, the perihelion of the moving body 
advances a fraction of a revolution equal to 

during each revolution. 

6. As regards the approximation to u, we write 

, _ _ TTCr 

(p'—'VJ <T, X — ~ — ;r~ j 

2pK 

and as k^ is small, the quantity (j of the elliptic functions is 

1 

Now from the expression for en -, we have 

on ~ = (1 — 4$sin^^')cos^, 
P 

so that we want cos.^^ But 



gm- 



— 1 -r O -r-r , 



and therefore 



m^' 



cos^ = cos 4 f 1*^3 — jo- 



.^, m^ 
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to the order retained, so that, 



a . f» w,^ . m^ 



en - •=: COS cr+3T^<7sin<7— --r€ sin^ cr cos o; 
p fi^ hf 

Again, from the expression for dn — , we have 

P 

dn — =» 1— Sgsin^a? 
p 

= 1 — 2 -r-r- € sin^ cr. 
^^ 

Consequently 

1 + g^ cr/p _ 1 1 "f cos cr + 3 m^/h^ or sin cr -> m^/ A^ e sin^ <r cos <r 
1 + ^^ o'/p 1 ■— m^/h^ e sin^ cr 

= i ( l + coscr + 3-~r crsincr + -T^^sin^a-) 
to our order of approximation. Also 

2 2 

so that ^ = g-(l_e); 

and 2m(y8~-7) = 2 -75 = 4 



a/B "" ^' 



so that ^ — ry = 2 ~- 6. 



Thus 



1^ = 7-„ ( 1 — e) + —^14- cos cr+ 3 -:r5 <x sm 0-4- — « ^ sm^ o* 
=== ~ ^ l + 6COScr + 3—:j-.crsmcr + -~g2gij^2^ 

— T2] l + gcos{^ — sy)4-T2^sin^(<^— t«y) 



AM '^ ' ' A2 



+ 3-p-(^-w)8in(^— w) k 



7. If a closer approximation is desired for any body in which the quantity 
m^/h^ is not of so small an order as in the known solar system, fuller values 
of the roots a, ^, 7 of U = should be found, and can easily be found. But 
an inspection of the value of 6^ + 4H^, upon which these roots depend, shows 
that, for any known bodies, the further approximations are ordinarily not 
worthy of account. 
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